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TRACK-TO-TRACK  ASSOCIATION  IN  OCEAN  SURVEILLANCE 


INTRODUCTION 

The  track-to-track  association  problem  in  ocean  surveillance  is  discussed  in  this 
report.  A family  of  procedures  is  described  for  deciding  if  two  observed  tracks  are 
from  the  same  target,  and  these  procedures  are  discussed  within  the  operational  context 
of  a large-scale  automated  surveillance-data  processing  system. 

The  Navy’s  roles  in  protecting  and  maintaining  the  sea  lines  of  communication  and 
in  supporting  other  areas  of  national  concern  require  a major  effort  in  the  surveillance  of 
ocean  surface  traffic.  This  effort  is  required  for  monitoring  possibly  hostile  combatant 
ships  to  protect  the  U.S.  fleet  and  also  for  such  activities  as  guaranteeing  the  freedom  of 
movement  of  U.S.-bound  oil  tankers  and  locating  and  warning  U.S.-flag  merchant  ships  of 
possible  dangers.  The  increasing  need  for  effective  and  timely  surveillance  information 
has  generated  a requirement  for  diverse  data  sources  and  for  the  capability  to  process 
surveillance  data  more  rapidly  than  can  be  done  by  unaided  analysts.  The  Navy  is  dev  sl- 
oping systems  which  accept  surveillance  data  of  various  types  from  a variety  of  sources, 
process  the  data  at  speeds  adequate  to  meet  increased  report  arrival  rates,  and  store  and 
maintain  the  finished  product  in  a readily  accessible  form.  A useful  format  for  stored 
information  is  a track,  defined  as  a time-ordered  sequence  of  positions,  with  all  sequence 
elements  relating  to  the  same  target.  Instances  of  the  track-to-track  association  problem 
occur  when  two  unidentified  tracks  are  examined  to  determine  whether  they  could  have 
been  developed  on  the  same  target  or  when  two  tracks  identified  as  relating  to  the  same 
target  are  examined  to  decide  if  the  identification  is,  in  both  cases,  correct. 

This  report  is  one  of  a series  of  studies  resulting  from  continuing  Naval  Research 
Laboratory  efforts  in  the  development  and  analysis  of  surveillance  systems  and  concepts. 
Recent  NRL  studies  of  surveillance  problems  include  Refs.  1 through  5.  For  the  case  of 
two  sensors,  each  of  which  has  observed  the  position  of  a vessel,  Ref.  5 developed  a 
method  for  determining  whether  one  target  or  two  targets  had  actually  been  detected. 

The  present  study  extends  this  situation  to  the  case  of  target  tracks.  The  test  procedures 
discussed  in  this  report  are  not  unique  to  the  track  association  problem.  Reference  6,  for 
example,  notes  that  similar  procedures  are  applicable  in  surveillance,  quality  control,  and 
motion  detection.  The  purpose  of  this  report  however  is  to  describe  the  relevance  of  each 
of  the  test  procedures’  operating  parameters  to  the  track  association  problem  and  to  indi- 
cate how  these  parameters  will  affect  operational  implementation  of  the  procedures. 
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RESULTS 

The  major  analytic  results  of  this  investigation  involve  the  resolution  of  track-to- 
track  association  problems  through  a series  of  statistical  tests.  The  basic  problem  of 
determining  whether  two  observed  tracks  relate  to  the  same  target  has  been  formulated 
as  a test  involving  the  chi-square  distribution.  An  extended  problem  of  determining 
which  of  two  candidate  tracks  a given  track  can  best  be  associated  with  has  been  formu- 
lated as  a test  of  the  chi-square  distribution  against  a noncentral  chi-square  alternative. 
Finally  the  extended  problem  with  the  added  feature  that  decisions  can  be  deferred 
pending  the  receipt  of  additional  data  has  been  formulated  as  a sequential  test.  The  test 
for  the  basic  problem  can  be  readily  implemented  within  an  automated  surveillance-data 
processing  system  and  can  be  carried  out  automatically  without  requiring  inputs  from  sur- 
veillance analysts.  The  other  tests  will  require  inputs  from  analysts.  The  estimates 
required  of  the  analyst  are  of  operationally  significant  parameters,  although  fortunately 
in  many  cases  the  final  decision  of  the  test  is  not  sensitive  to  moderate  changes  in  the 
parameters. 


OPERATIONAL  SITUATION 

For  this  study,  processing  of  surveillance  information  consists  of  maintaining  and 
updating  data  files  on  tracked  targets.  The  procedure  is  assumed  to  be  carried  out  by  a 
processing  system  which  .includes  surveillance  analysts,  computers,  and  analyst  aids,  both 
manual  and  computer-driven.  Computers  perform  the  bulk  of  the  routine  data-processing 
functions,  freeing  the  analysts  for  making  decisions  in  difficult  cases  and  for  providing 
judgments  which  are  beyond  the  capabilities  of  automated  processes. 


Track  association  schemes  will  be  used  to  associate  tracks  resulting  from  a report-to- 
track  association  process.  These  tracks  must  contain  the  best  available  information  on 
observed  target  movements.  Report-to-track  association  processes  operate  on  sets  of  re- 
ported target  positions  and  sets  of  established  tracks  and  attempt  to  match  correctly  the 
elements  of  these  sets.  The  pairing  of  a reported  position  with  an  established  track  ex- 
tends the  track  to  a new  position,  either  the  position  observed  in  the  report  or  a smoothed 
position  obtained  from  operations  on  the  reported  position  and  track  parameters.  When 
two  or  more  tracks  exhibit  similar  positions  over  an  interval  of  time,  a track  association 
scheme  is  called  into  use. 


j 
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The  tracks  upon  which  a track  association  scheme  is  to  be  used  will  contain  the 
underlying  positions  and  reported  times,  together  with  the  covariance  matrices  associated 
with  these  positions.  If  the  track  comprises  independently  observed  positions,  then  the 
covariance  matrix  will  indicate  no  correlation  between  position  coordinates  for  different 
observation  times.  However,  if  the  track  is  made  up  of  smoothed  positions,  then  nonzero 
correlations  will  occur.  This  is  because  track  smoothing  schemes  are  based  on  estimating 
target  positions  by  applying  a computational  procedure  either  to  previous  smoothed  posi- 
tions or  directly  to  previous  observed  positions,  so  that  each  smoothed  position  contains 
information  about  previous  smoothed  positions;  thus  the  covariance  matrix  associated 
with  a set  of  smoothed  positions  will  contain  entries  indicating  correlation  between  the 
position  estimates. 
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In  this  report  the  surveillance  information  is  considered  to  be  generated  by  two  in- 
dependently operating  surveillance  systems.  The  first  system,  called  the  prim  -try  system, 
is  defined  as  possessing  superior  capabilities  in  a number  of  important  system  haracteris- 
tics,  such  as  high  detection  probability,  high  frequency  of  observation  of  a giver,  area,  or 
small  position  estimate  errors.  The  other  system  (the  secondary  system)  is  composed  of 
subsystems,  no  one  of  which  is  as  effective  as  the  primary  system.  The  members  of  the 
secondary  system  may  have  poor  detection  capabilities,  may  detect  targets  only  through 
infrequently  exposed  characteristics,  or  may  possess  large  localization  errors.  They  may 
however  possess  important  characteristics  which  the  primary  system  does  not.  For  exam- 
ple the  secondary  system  might  be  able  to  identify  the  target  by  nationality,  by  class 
(merchant  or  combatant),  or  by  name.  Although  tracking  unknown  targets  with  the 
primary  system  may  provide  valuable  information  on  target  movements,  the  addition  of 
information  on  target  identity  from  a secondary  source  will  increase  the  information 
content  of  both  tracks  and  may  provide  essential  data  for  analyzing  tactics  or  predicting 
intentions. 

It  is  assumed  that  throughout  the  major  part  of  the  surveillance  operation  the  data 
from  the  two  surveillance  systems  are  processed  independently  and  the  report-to-track 
association  and  track  continuation  processes  that  were  described  are  conducted  separately 
for  the  two  systems.  Thus  prior  to  the  employment  of  the  track-to-track  association 
process  the  data  base  of  the  surveillance  processing  system  contains  tracks  generated  by 
the  primary  system  and  tracks  generated  by  the  secondary  system.  At  certain  times 
however  the  track  families  are  examined  to  determine  whether  any  of  them  refer  to  the 
same  target  and  thus  can  be  merged  into  a single  track. 


BASIC  MODEL 

In  this  section  the  basic  model  employed  in  the  investigation  is  described;  this  model 
involves  one  track  from  the  primary  system  and  one  track  from  the  secondary  system. 

The  primary  system  is  denoted  , and  its  track  is  denoted  T1 ; the  secondary  system  is 
denoted  S0,  and  its  track  is  denoted  T0. 

Figure  1,  which  shows  Tq  and  Tj,  is  representative  of  the  basic  problem  formula- 
tion. The  tracks  have  been  generated  by  connecting  target  positions  as  generated  by  the 
two  sensor  systems.  In  the  basic  model  each  track  consists  of  p positions,  and  the  two 
systems  are  assumed  to  have  both  made  their  observations  at  the  same  times  t\>  t2,  •••, 
tp.  The  ith  position  of  fraek  Tu  is  denoted  by  the  coordinate  pair  (x\h\  y\!^),  for 
i = 1,  2,  ...,  p and  k * 0,  1.  Thus  each  can  be  associated  with  the  n-by-one  (where 
n = 2p)  column  vector 
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It  is  assumed  that  the  n-tuples 


x(k ) V(A) 
\X1  ’ yl  > 


«’) 


are  observations  from  an  n-dimensional  normal  distribution.  For  an  arbitrary  random 
variable  x,  let  x represent  the  expectation  of  x.  Then  each  observed  track  is  regarded 
as  a sample  from  an  n = 2p  dimensional  normal  distribution 


where  each  is  a 2p-by-2p  covariance  matrix.  The  elements  of  the  matrices  2^  are 
determined  by  the  characteristics  of  the  sensor  systems  and  possibly  by  the  analytic 
methods  used  to  obtain  smoothed  coordinates  (x{k\  yj^))  from  the  reported  observed 
positions. 


The  test  statistic  to  be  used,  denoted  I?2,  is  defined  by 


(1) 


where  the  prime  indicates  matrix  transpose.  Given  2q  and  Sj,  it  is  known  [7]  that 
there  exists  a nonsingular  matrix  P such  that  P(T.q  + '£\)  P'  - I,  where  I is  the  ti-by-n 
identity  matrix.  It  will  be  shown  that  the  statistic 
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has  a distribution  whose  characteristics  are  useful  in  this  investigation.  Since  it  can  be 
showr<  that  W2  = R2,  in  the  following  we  will  consider  that  R 2 = \P(Tq  - Tx ) I2  and  use 
the  notation  r2  for  observed  values  of  R2  or  for  values  which  have  been  calculated  from 
a specified  set  of  track  positions. 


NULL  HYPOTHESIS  H0 

The  null  hypothesis  to  be  tested  is  H0 : T0  and  Tx  are  from  the  same  target.  More 
specifically,  let  represent  the  column  vector  of  expected  values 


Then  an  equivalent  formulation  of  the  null  hypothesis  is 

#o:  70  = T'x 


~ = A 

where  0 is  the  2p-by-l  column  vector  with  all  components  equal  to  0. 


*8 


Tim* 

Fig.  1— The  track  associatt  n problem  in  the  basic  model,  in  which  the  primary  system 
(subscript  1)  and  secondary  system  (subscript  0)  observe  track  positions  simultaneously 

5 


mciiWjfTM  irr-i 


wr 


7rr  *5r  snvmw. 


gjgyyf  yyww  yffly 


HOWARD  L.  WIENER 


Since  the  tracks  Tq  and  T1  are  samples  from  n-dimensional  normal  distributions, 
with  n = 2p,  the  difference  T0  - Tj  is  a sample  from  an  n-dimensional  normal  distribu- 
tion with  expectation  T0  - 7\  and  covariance  matrix  2 = 2q  + Sj.  As  mentioned,  under 
H0,  Tq  - Ti  = £.  With  the  matrix  P defined  as  just  before  equation  (2),  it  follows  by 
properties  of  the  multivariate  normal  distribution  [8,  pp.  24-25]  that  P(T0  - Tx ) is  a sam- 
ple from  an  n-dimensional  normal  distribution  with  zero  mean  and  identity  covariance 
matrix.  Thus  [7,  p.  56]  under  H0,  R2  has  the  (central)  chi-square  distribution  with 
n = 2p  degrees  of  freedom,  with  probability  density  function  (p.d.f.) 

f0(r2;n)  = [2Pr(p)]P~1  e^12.  (3) 

The  null  hypothesis  H0  is  rejected  whenever  the  computed  value  r2  of  R2  seems  too 
large.  A type-I  error  occurs  when  the  null  hypothesis  is  incorrectly  rejected;  the  proba- 
bility of  this  error,  the  level  of  significance,  is  denoted  a.  For  a specified  value  of  a and 
a given  number  of  observed  positions,  the  acceptance  region  for  H0  is  bounded  below  by 
0 and  above  by  the  value  x2  satisfying  the  expression 


Values  of  x2  for  given  a and  p are  found  in  tables  of  the  chi-square  distribution,  such  as 
in  Ref.  9;  Table  1 contains  representative  values.  For  computer  application  either  the 
chi-square  tables  can  be  stored  in  memory  or  analytic  expressions  can  be  used  to  compute 
approximations  to  x2;  such  approximations  are  discussed  in  Chapter  17  of  Ref.  10  (Vol.  1). 

It  is  valuable  to  consider  the  situation  that  inspires  the  rejection  of  the  null  hypoth- 
esis Hq  ; these  considerations  lead  to  discussions  of  when  the  significance  test  ought  to  be 
carried  out  on  two  observed  tracks.  It  is  clear  that  the  test  should  be  carried  out  when 
and  only  when  there  is  doubt  as  to  the  source  of  two  obse„"ed  tracks.  As  an  example  of 
a situation  in  which  the  test  is  not  required,  suppose  there  are  tracks  T0  and  T1  as  shown 
in  Fig.  1,  where  the  primary  surveillance  system  has  a high  detection  probability.  For  II  q 
to  be  rejected,  the  track  T0  must  be  assumed  to  have  come  from  a target  other  than  that 
generating  I\ . However,  as  shown  in  the  figure,  no  such  alternate  track  has  been  detected, 
and  the  existence  of  an  alternate  target  is  a low  probability  event  in  the  light  of  the  pri- 
mary system’s  high  detection  probability.  Thus,  since  no  other  target  could  have  gener- 
ated Tq,  it  is  sufficient  to  associate  the  two  tracks  without  recourse  to  the  test. 

It  is  not  sensible  to  carry  out  the  test  unless  rejection  of  HQ  is  operationally  feasible. 

If  the  situation  at  hand  does  not  admit  the  possibility  of  at  least  one  additional  target  for 
association  with  Tq,  then  the  two  candidate  tracks  should  be  associated.  Situations  in 
which  the  test  should  be  carried  out  will  always  admit  the  possibility  of  separate  targets 
underlying  the  tracks  T0  and  Tl . For  example,  whenever  the  primary  surveillance  sys- 
tem has  only  a moderate  detection  probability,  it  is  possible  that  a target  may  be  present 
yet  not  be  detected  sufficiently  often  to  have  generated  a track.  Another  possibility  is 
that  a number  of  candidate  tracks  have  been  detected  by  the  primary  system,  so  that 
should  the  test  fail  there  will  be  other  candidates  for  association  with  T0.  This  case  is 
discussed  in  the  following  section. 
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Table  1— Acceptance  Regions  for  the  Hypothesis  HQ 


Number  p of  Observed 
Positions  in  Each  Track 

Maximum  Value  x2  for  Acceptance  of  HQ 

a = 0.30 

a = 0.20 

a = 0.10 

a = 0.05 

a = 0.01 

2 

4.878 

5.989 

7.779 

9.488 

13.277 

4 

9.524 

11.030 

13.362 

15.507 

20.090 

6 

14.011 

15.812 

18.549 

21.028 

26.217 

8 

18.418 

20.465 

23.542 

26.296 

32.000 

10 

22.775 

25.038 

28.412 

31.410 

37.566 

15 

33.530 

36.250 

40.256 

43.773 

50.892 

ALTERNATE  HYPOTHESIS  Hx 

Although  the  track  association  test  should  be  performed  only  when  it  is  reasonable 
to  suggest  that  track  T0  could  be  related  to  a target  other  than  the  source  of  Tx , it  is 
not  always  the  case  that  a definite  alternate  source  for  T0  can  be  proposed.  It  is  one 
thing  to  accept  the  possibility  of  an  alternate  source,  but  it  is  another  to  specify  one. 

In  the  situation  discussed  in  this  section,  an  observed  track  is  specified  as  relating  to  a 
possible  alternate  source  for  TQ.  In  this  case  one  can  specify  an  alternate  hypothesis  Hx 
and  determine  the  power  of  the  track  association  test  against  this  alternative. 

Figure  2,  adapted  from  Ref.  11,  illustrates  an  operational  situation  in  which  this 
case  could  arise.  This  figure  contains  examples  of  ship  track  histories  a surveillance 
system  such  as  the  primary  system  might  obtain  over  an  interval  of  time  for  a given  area 
of  interest.  Each  of  these  tracks  represents  a separate  target.  This  example  contains 
both  ships  which  appear  to  be  moving  randomly  and  ships  which  move  in  air  apparently 
more  predictable  manner.  The  illustration  shows  much  track  crossing  and  few  cases  in 
which  two  or  more  tracks  remain  pa.  .rllel  for  an  extended  period.  Suppose  that  a track 
T0  from  the  secondary  system  appears  to  be  related  to  one  of  the  tracks  in  this  illustra- 
tion. If  the  association  test  should  fail,  then  it  would  appear  that  there  will  be  at  most 
one  other  track  with  which  T0  might  be  associated.  If  the  primary  system  has  a relatively 
high  detection  probability,  then  it  is  unlikely  that  it  would  not  have  detected  the  target 
which  generated  T0.  Consequently  it  is  assumed  that  there  will  be  another  candidate 
track  for  T0  and  thus  that  if  association  with  the  first  track  should  fa;l,  then  association 
with  the  second  track  should  be  accepted. 

Figure  3 illustrates  the  extension  of  the  basic  model  to  the  case  of  testing  the  hy- 
pothesis Hq  against  a specific  alternative.  The  tracks  T0  and  Tx  are  as  before.  Another 
track  T2  is  now  present  and  is  assumed  to  have  been  obtained  from  the  primary  system 
Sx . The  track  TQ  must  relate  to  one  of  the  other  tracks;  if  the  hypothesis  that  T0  is 
associated  with  Tx  is  rejected,  then  T0  will  be  associated  with  T2.  The  statistical  distribu- 
tions associated  with  T0  and  Tx  are  as  was  described  in  the  section  on  the  basic  model. 
The  observations  comprising  T2  constitute  a column  vector 
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and  this  vector  is  assumed  to  be  a sample  from  a 2p-dimensional  normal  distribution  with 
expectation 


and  2p-by-2p  covariance  matrix  22.  It  is  clear  that  this  case  admits  a specified  alternate 
hypothesis  Hi , which  states  that  T0  is  related  to  T2 . The  alternate  is  given  by 

Hl:  T0  = T2  . 


Distribution  of  R 2 Under  H^;  Power  of  the  Test 


Under  Hi  the  difference_vector  T0  - is  a sample  from  a 2p-dimensional  normal 
distribution  with  expectation  T2  - and  covariance  matrix  £ = E0  + 2), . Thus  [7,  p. 
56]  R 2 has  the  noncentral  chi-square  distribution  with  n = 2p  degrees  of  freedom  and 
noncentrality  parameter 


mmmmm 
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where  matrix  P is  as  defined  previously.  The  probability  density  frunction  of  the  noncen- 
tral chi-square  distribution  with  n degrees  of  freedom  and  noncentraiity  X is  defined  by 


1 /r2 

Mr2;*,  X)  = ~(j 


(»-2)/4  ^_<r2+X)/2 


I(«-2)/2^rV^»  r > 0, 


where  1% (...)  is  the  modified  Bessel  function  of  the  first  kind  and  order  K [Ref.  10].  Thus 
for  a given  set  of  p positions,  a specified  value  for  the  level  of  significance  a,  and  specified 
value  of  X,  the  probability  of  incorrectly  accepting  the  hypothesis  Hq  when  is  true  is 

X2 

13=  f M*2;".  X)dr2  , (5) 

Jo 

where  n = 2p  and  where  the  upper  limit  x2  is  defined  by  expression  (4).  The  power  of 
the  test,  n = 1 - 0,  is  given  by 


n “ J fl  (r2 ; n,  X)  dr2  . (6) 

X2 

For  track-to-traCk  association  the  test  of  H0  against  H1  is  a test  of  the  central  chi- 
square  distribution  against  a noncentral  chi-square  alternative.  Tables  relating  to  this  test 
are  available  [for  example,  Ref.  12].  Reference  13  describes  a rapid,  compact  program 
for  carrying  out  computations  based  on  the  noncentral  chi-square  distribution.  This  pro- 
gram is  suitable  for  implementation  within  automated  tracking  algorithms  or  for  use  within 
computer-driven  analyst  aids. 

Figure  4 indicates  the  power  of  the  test  for  levels  of  significance  of  0.10  and  0.01 
and  for  situations  involving  six  and  15  observed  positions  in  each  track.  When  all  other 
parameters  are  held  fixed,  the  power  of  the  test  increases  if  X increases,  increases  if  the 
specified  level  of  significance  increases,  and  decreases  if  the  number  of  observed  positions 
increases.  Generally  speaking,  the  operational  situation  will  determine  the  values  of  p, 

Tq,  Ti  T<i , Z0>  and  The  error  probabilities  a and  (3  can  be  specified.  However, 
the  value  of_X  must  be  estimated  by  the  surveillance  analyst,  since  the  value  of  the  differ- 
ence vector  T2  - Tj  cannot  be  assumed  to  be  known.  An  important  question  underlying 
the  test  procedure  is  the  sensitivity  of  the  resulting  errors  a and  (3  to  variations  in  the  X 
value  employed.  It  must  be  determined  whether  the  test  procedure  requires  great  accu- 
racy on  the  part  of  the  analyst. 


) 

i 


Sensitivity  to  Estimates  of  X 

liiv  oarameter  X is  determined  by  the  difference  vector  T2  - T ^ and  by  the  matrix 
sum  Z = r0  + 21;  since  the  latter  factor  will  generally  be  known,  the  requirement  that 
X be  known  reduces  to  a requirement  that  the  difference  vector  be  estimated,  since  it 
cannot  be  known  exactly.  The  effects  of  variation  in  either  the  difference  vector  or  in 
X on  the  results  of  the  test  procedure  can  be  ascertained  only  by  a detailed  parametric 
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analysis  of  the  relations  between  p,  a,  j3,  T2  - T1(  and  2.  Some  insights  into  the  general 
problem  of  test  sensitivity  to  X can  be  obtained  from  an  analysis  of  the  two-candidate 
test  in  the  following  special  case.  In  the  two-candidate  track  association  problem  the 
choice  of  which  track  to  use  as  a basis  for  the  null  hypothesis  Hq  is  arbitrary.  Therefore 
in  the  absence  of  any  supplementary  information  it  will  be  assumed  that  errors  in  failing 
to  make  a correct  association  are  as  important  for  the  alternate  hypothesis  as  for  the  null 
hypothesis.  In  this  case  it  is  reasonable  to  set  a = (3.  Having  specified  a,  one  can  proceed 
as  before  by  finding  that  value  of  x2  such  that  equation  (4)  holds.  Using  this  value  as 
the  test  criterion,  one  can  then  employ  equation  (5)  to  determine  that  value  of  X such 
that  the  resulting  (3  value  is  equal  to  ot.  Analysis  of  the  resulting  relations  between  a(  =/3) 
and  X can  indicate  how  precisely  X must  be  estimated  to  obtain  a test  with  level  of  signif- 
icance a and  power  1 - a. 

Figure  5 shows  the  sensitivity  of  the  resulting  test  errors  to  estimates  of  X by  ex- 
hibiting contours  along  which  a = (3  for  p values  of  6,  15,  and  24.  Since  the  power  of 
tire  test  is  an  increasing  function  of  X,  this  figure  shows  that  to  guarantee  a specified 
maximum  error  probability  it  is  necessary  only  to  guarantee  that  the  actual  value  of  X 
exceeds  a stated  threshold.  For  example,  for  the  case  of  15  observed  track  positions, 
if  one  can  guarantee  that  X is  greater  than  30,  then  the  test  will  achieve  error  rates  of 
less  than  8 percent.  The  lower  threshold  for  estimates  of  X increases  as  the  number  of 
observed  positions  increases  and  as  the  allowable  maximum  error  probability  decreases. 
Therefore,  although  all  that  the  procedure  requires  is  a decision  that  X exceed  some 
stated  minimum  value,  it  may  be  difficult  to  assure  that  this  decision  is  correct. 


Fig.  4— Power  of  the  test  of  Hq  against  H1. 
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NONCENTRALITY,  X. 


Fig.  !>— Sensitivity  of  the  test  of  HQ  against  to  changes  in  the 
noncentrality,  for  a = (3. 


An  example  will  illustrate  the  conceits  involved  in  relating  requirements  on  X to 
requirements  on  the  estimated  value  of  T2  - T1 . Suppose  that  tracks  Tx  and  T2  were 
observed  with  a sensor  which  generated  a covariance  matrix  with  no  off-diagonal  elements 
and  such  that  both  a2  and  a2  equaled  1 n.mi.2  for  each  observed  position.  Suppose  that 
T o was  generated  by  a system  whose  covariance  matrix  also  was  diagonal,  but  with 
a2  = Oy  = 4 n.mi.2  for  each  position.  Finally,  suppose  six  positions  (12  coordinates) 
have  been  observed  and  that  an  analyst  can  estimate  that  each  element  in  the  difference 
vector  f2  - T1  is  at  least  2.0  n.mi.  Then  the  minimum  value  for  X will  be  9.6,  and  one 
can  conclude  from  Fig.  5 that  the  test  could  be  subject  to  errors  greater  than  20  percent. 
Suppose  on  the  other  hand  that  errors  of  at  most  5 percent  were  required  for  tracks 
comprising  15  observed  positions.  This  would  require  a value  of  X of  at  least  35,  and 
with  the  covariance  matrix  structure  described  it  would  be  necessary  to  be  correct  in 
deciding  that  the  elements  in  the  difference  vector  were,  on  the  average,  at  least  2.4  n.mi. 


Computer  Implementation 

Computer  implementation  of  a procedure  based  on  these  results  could  be  readily 
achieved  as  an  interactive  analyst  aid.  Since  the  majority  of  the  required  parameter 
values  would  already  be  available,  an  analyst  would  have  to  enter  his  best  judgment  of 
the  magniture  of  the  elements  in  the  underlying  mean  difference  vector.  The  program 
would  then  generate  the  test  procedure,  identifying  the  operable  value  of  a,  computing 
tire  value  of  the  upper  limit  x2.  computing  the  value  of  r2,  and  determining  whether  or 
not  r2  exceeded  x2>  It  would  respond  to  the  analyst  with  a statement  of  the  decision 
which  was  made,  together  with  an  indication  of  the  error  probabilities  present  in  the  test 
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i used.  Conversely,  if  an  analyst  were  to  specify  a desired  level  of  significance  a,  the  pro- 

i gram  could  respond  with  the  minimum  value  of  X required  to  guarantee  an  error  P no 

greater  than  a. 


1 

! 

I 

i- 

j| 

*■ 

\ 

ij 

*T. 

1. 


FURTHER  ASPECTS  OF  THE  PROBLEM 

In  tliis  section  the  preceeding  discussion  is  extended  to  three  additional  topics.  The 
first  is  a method  for  treating  the  situation  when  the  basic  assumption  of  simultaneous 
observations  is  relaxed.  The  second  is  the  case  in  which  immediate  association  decisions 
are  not  always  required  and  can  be  deferred  pending  the  receipt  of  further  data.  Finally 
the  third  is  the  importance  of  selecting  proper  limits  for  the  error  probabilities  a and  p 
in  the  tests  and  the  consequences  of  making  track  association  errors  within  a surveillance- 
data  processing  system. 


Nonsimultaneous  Observations 

The  basic  track-to-track  association  model  assumed  that  the  positions  constituting 
the  two  tracks  were  observed  at  the  same  times.  Usually  different  surveillance  systems 
will  make  detections  at  different  times;  hence  for  any  pair  of  generated  tracks  a set  of 
nonsimultaneous  observations  is  more  likely  than  a set  of  simultaneous  observations.  The 
following  discussion  is  a description  of  a method  of  testing  tracks  with  nonsimultaneous 
observations  for  equality  of  means.  The  method  transforms  the  case  of  nonsimultaneous 
observations  to  the  case  of  simultaneous  observations.  This  transformation  involves  the 
generation  of  interpolated  points  within  each  track  so  that  both  of  the  resulting  new 
tracks  will  have  contituent  points  associated  with  the  same  set  of  observation  times.  The 
set  of  observation  times  contains  all  times  at  which  either  of  the  two  original  tracks  was 
observed. 

Figures  6 and  7 represent  the  situation.  Figure  6 shows  track  T1 , observed  at  times 
t2,  f3,  f4,  and  t6,  and  track  T0,  observed  at  times  fly  t5,  and  t7.  To  test  the  hypothesis 
that  these  tracks  relate  to  the  same  target,  two  new  tracks  Tf,  and  T*  have  been  created 
(Fig.  7),  both  of  which  comprise  all  the  observation  times  for  both  Tx  and  T0.  The 
points  in  the  new  tracks  are  of  two  types.  First,  all  points  originally  observed  for  a track 
are  contained  in  the  related  new  track.  Second,  if  an  original  track  Tk  contains  two  se- 
quential points,  observed  at  times  tj  and  tj,  and  if  the  other  track  contains  a point  observed 
at  an  intermediate  time  t,  where  t;  < t < tj,  then  an  intermediate  point  is  interpolated 
between  the  original  points.  Letting  D = (t  - - ff),  the  new  point  is  assumed  to 

occur  at  time  t and  to  be  located  at  a fraction  D of  the  way  between  the  originally 
observed  points.  If  xt  and  Xj  are  the  jc  coordinates  of  Tk  at  times  t,-  and  tj  respectively 
and  x*  is  the  jc  coordinate  of  the  new  interpolated  point  in  track  T%  at  time  f,  then 

x*  = Xj  + D(xj  - xt) 

_ (tj  - t)  (t  - tj) 

(tj  - tj)  Xi  (tj  - tj)  Xj- 
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Fig.  7— Nonsimultaneous  tracks  of  Fig.  6 after  transformation  to  the  case  of 
simultaneous  observations 


These  three  points,  the  two  original  ones  and  the  new  interpolated  one,  are  all  made 
members  of  the  new  track  Tfy  The  process  is  carried  out  for  both  of  the  original  tracks 
and  all  of  the  time  points  for  which  the  interpolation  is  feasible,  resulting  in  the  situation 
shown  in  Fig.  7.  In  this  figure  the  two  new  tracks  T-f  and  Tq  both  have  “observation 
points”  at  all  times  f2  through  t6.  Since  tires e points  have  been  obtained  by  a linear 
transformation  of  normal  random  variables,  their  related  covariance  structures  can  be  ob- 
tained [8]  and  they  may  be  used  as  data  for  the  test  of  the  null  hypothesis  H0  described 
for  the  basic  model.  The  one-point  extensions  in  track  Tq,  both  at  the  start  and  at  the 
end,  are  not  used  in  the  hypothesis-testing  computations,  as  the  interpolation  scheme 
does  not  generate  extensions  from  track  T1 . 

This  extension  of  the  basic  model  to  the  case  of  nonsimultaneous  observations  has 
been  programmed  for  use  within  an  automated  track  correlation  process.  This  program 
is  described  in  Appendix  A and  listed  in  Appendix  B. 


Sequential  Tests 

A situation  was  discussed  in  which  track  T0  was  to  associate  with  one  of  two  can- 
didate tracks  7\  and  To  as  a result  of  a testing  procedure.  It  is  not  always  the  case  that 
an  immediate  track  association  decision  is  required;  quite  often  the  problems  inherent  in 
making  association  errors  are  severe  enough  to  require  that  association  decisions  not  be 
made  without  sufficient  supporting  evidence.  One  method  of  accomplishing  this  is 
through  the  use  of  sequential  tests  based  on  reexamination  of  the  data  set  each  time  new 
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data  elements  are  obtained.  In  the  following  discussion  the  structure  of  sequential  testing 
procedures  within  the  context  of  the  track  association  problem  is  described  and  some 
problems  inherent  in  implementing  the  procedure  are  indicated. 

As  described  in  Refs.  14  and  15,  sequential  testing  procedures  are  useful  in  deciding 
whether  observed  data  correspond  to  a null  hypothesis  Jfi0,  whether  they  correspond  to  a 
specified  alternate  hypothesis  H1 , or  whether  additional  data  should  be  obtained  before 
making  a final  decision.  Let  g0  n and  g\  n be  the  probability  density  functions  of  the 
multivariate  data  set  (jcj,  x2,  ...,  xn)  at  the  nth  step  of  the  process,  under  hypothesis  Hq 
and  Hi  respectively,  let 

^ x2'  “•>  xn ) 

gO  ,n(*l»  *2  * ) 

and  suppose  that  the  error  probabilities  a and  0 have  been  specified.  The  sequential  test- 
ing procedure  is  then  defined  by  the  following  criteria: 


' i 

, 1 

, i 

lthn  < 1 _a’ 

accept  Hq  ; 

(7a) 

,s 

.t> 

i 

'!  1 
» 

! 

i 

if  < hn  < ^ , 

l-o  n a 

obtain  another  observation 

(7b) 

\ 

: $ 

\ • 

i 

n 

if^—^  < hn  , 
a 

accept  Hi  . 

(7c) 

\.,fu 

Note  that  the  criteria  do  not  depend  on  the  sample  size  n.  These  criteria  were  based  on 
the  assumptions  that  the  successive  observations  were  stochastically  independent  samples 
and  that  the  sequential  procedure  will,  with  probability  1,  eventually  terminate.  For  a 
rich  family  of  situations  the  assumption  of  stochastically  independent  samples  can  be  re- 
laxed and  the  procedure  will  still  be  valid,  it  has  been  assumed  here  that  the  conditions 
underlying  the  track  association  problem  are  such  that  the  sequential  test  procedure  is 
valid. 

For  track  association  the  assumptions  of  underlying  normal  distributions  imply  that 
under  the  null  hypothesis  H0  specified  earlier  T0  - Tl  has  tire  density  function 

&Q ,n  ~ So(To  -Tx\p) 

’ (5ijpTiTT«  cxp  [-i^-n^ro-TO 

and  under  H1  specified  earlier  T0  - Tj  has  the  density  function 


1 imuiirmruanmi 


15 


HOWARD  L.  WIENER 


is 

l 

'] 

{ 

i 

t 


l 

) 

i 


i 


i 


8l,n  = 8x(T0-Ti;P,  T2-Tx) 


{2n)P  IS |i/2  exP^2{[(T°-Tl)  - (Ta-TOl'Z-HCTo-Ti)  - (T2-Tx)]} 
It  follows  that 


__  8\,n 
80/> 


= exp 


--HTo-TO'^-^Ta-ri)  - (Tz-T^'Z-HTo-T!) 


+ (T2  - T\)'  £-1(T2  -Ti)] 


= exp 


jj[X-2(T0-T1)'I-1(T2-:T1)]J. 


As  discussed  previously,  it  is  reasonable  to  consider  that  the  error  probabilities  a and  ft 
are  equal.  Thus,  with  the  notation 

2 = X - 2(r0-T1)'£-1(T2-T1), 

criteria  (7)  reduce  the  foil  <wmg: 

if^r^F<z<te(ir)2- 


accept  Hq  , 

(8a) 

wait  for  additional  data, 

(8b) 

accept  hi  . 

(8c) 

The  performance  of  this  procedure  is  determined  by  the  relations  between  tire  ob- 
served values  T0  - T1  and  2 and  the  estimated  mean  difference  vector  T2  - 7’, . Detailed 
investigations  based  on  realistic  predictions  of  the  observed  parameter  values  will  be  re- 
quired to  ascertain  the  sensitivity  of  the  sequential  procedure  to  estimates  of  the  differ- 
ence ve-  tor  in  an  operational  setting.  At  this  point  however  examination  of  the  criteria 
(8)  permit  general  conclusions  to  be  made  regarding  the  test’s  performance.  Figure  8 
illustrates  them  criteria.  If  the  permissible  error  probabilities  are  low,  the  hypothesis 
acceptance  regions  are  reduced;  the  sequential  procedure  will  most  frequently  decide  to 
wait  for  additional  data  unless  a value  of  Z is  obtained  which  is  large  in  absolute  value. 

If  the  test  is  more  lenient— reflected  in  higher  allowable  error  probabilities— the  sequential 
procedure  will  more  frequently  result  in  a hypothesis  acceptance,  unless  the  absolute 
value  of  Z is  small. 
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ERROR  PROBABILITIES  a « (3 
FiH-  8—  Sequential  test  criteria  (8)  for  track  association. 


In  summary,  the  use  of  sequential  testing  procedures  may  be  valuable  when  two 
tracks  7\  and  7’2  are  candidates  for  association.  This  test  permits  additional  information 
to  be  collected  if  the  available  information  is  adequate  for  a hypothesis  acceptance  deci- 
sion. The  test  procedure  requires  the  estimation  of  the  mean  difference  vector  T2  - 711 . 
The  sensitivity  of  the  test’s  performance  to  estimates  of  this  value  is  i function  of  the 
observed  data  T0  - T1  and  covariance  matrix  X ; further  detailed  investigations  of  this 
factor  will  be  required  to  determine  the  worth  and  feasibility  of  implementing  the  sequen- 
tial procedure  within  an  operational  surveillance  processing  system. 


Specifying  the  Allowable  Error  Probabilities 

The  tests  described  in  this  report  require  the  specification  of  maximum  allowable 
probabilities  a and  p of  incorrect  decisions.  These  errors  have  real-world  costs,  h hose 
responsible  for  implementing  these  tests  in  an  operational  surveillance-data  processing 
system  must  analyze  the  consequences  of  making  incorrect  decisions  in  order  to  specify 
a and  fi. 

In  the  case  of  the  basic  model  the  only  decision  is  whether  or  not  to  associate  a 
given  pair  of  tracks;  as  developed,  the  test  requires  the  specification  of  a,  the  probability 
of  incorrectly  rejecting  the  proposed  association  of  the  two  tracks.  Errors  of  this  type 
will  require  the  surveillance  processing  system  to  maintain  multiple  tracks  in  situations 
when  maintenance  of  single  tracks  would  suffice.  The  data  base  will  thus  become  clut- 
tered with  redundant  tracks,  and  any  computational  routines  that  are  performed  on  each 
track  in  the  data  base  will  be  performed  more  than  necessary.  More  important  is  the 
possibility  of  losing  information  on  special-interest  targets  which  correct  track  association 
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would  have  provided.  For  example,  without  si’  porting  information  on  target  identities, 
position-only  tracks  must  all  be  subject  to  the  same  degree  of  attention  by  analysts.  How- 
ever, if  position-only  tracks  can  be  associated  with  trac!;s  carrying  identification  informa- 
tion, then  known  combatant  tracks  can  be  specified  for  special  attention  and  tracks  from 
noncombatants  can  be  specified  for  lesser  attention.  Thus  a low  false  rejection  rate  for 
Hq  will  increase  data-processing  efficiency  and  increase  the  value  of  the  information  in 
the  stored  tracks.  Still  the  consequences  of  falsely  associating  unrelated  tracks  can  be 
severe.  Whereas  false  rejection  results  in  the  loss  of  useful  information,  false  acceptance 
surely  results  in  the  generation  of  misleading  information.  Misclassification  of  a combat- 
ant as  a merchant,  or  conversely,  can  certainly  create  false  pictures  of  tactical  situations. 
Track  files  containing  incorrectly  associated  tracks  can  easily  serve  as  bases  for  future 
compounded  errors,  thus  leading  to  a possibly  useless  data  base.  Extensive  lists  of  the 
consequences  of  making  eithe?  i,  >c  of  error  can  be  generated  indefinitely,  but  these  ex- 
amples should  suffice  to  indier-e  ine  importance  of  proper  assessment  of  the  effects  of 
the  errors  on  the  operation  and  ultimate  value  of  the  data,  processing  system. 


SUMMARY  REMARKS 

This  report  has  discussed  a series  of  statistical  tests  for  attacking  the  track-to-track 
association  problem  in  ocean  surveillance.  The  tests  are  straightforward  and  can  be 
readily  implemented  within  many  automated  surveillance-data  processing  systems.  Com- 
putational algorithms-.exist  whicn  permit  rapid  computation  of  any  of  the  test  statistics 
involved;  consequently  these  test  procedures  can  he  implemented  “on  line,”  within  an 
automatic  processing  system,  or  they  can  bo  used  within  interactive  programs  to  provide 
quick  response  to  analyst  queries.  The  demand  for  inputs  from  analysts  will  be  at  a low 
level;  optimally,  analysts  may  be  asked  simply  to  judge  whether  a given  parameter  is  or  is 
not  within  a specified  range.  A major  value  of  these  tests  lies  in  their  permitting  greater 
use  of  the  various  types  of  information  contained  in  different  tracks  of  the  same  target. 
The  result  of  successful  implementation  of  the  tests  is  a data  base  which  is  low  in  errors 
and  low  in  redundancy. 
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Appendix  A 

A TRACK-TO-TRACk  ASSOCIATION  PROGRAM 


INTRODUCTION 

A program  for  track-to-track  association  tes's  is  described  in  this  appendix  and  listed 
in  Appendix  B.  The  program  is  designed  to  operate  on  tracks  which  are  defined  by  sets 
of  positions  and  which  are  accompanied  by  the  underlying  covariance  matrices.  The 
track  positions  may  comprise  locations  reported  by  a surveillance  sensor  or  may  consist 
of  estimated  positions  obtained  by  applying  a smoothing  scheme  Jo  a set  of  observations. 
In  the  former  case  the  covariance  matax  would  probably  be  block  diagonal,  with  two-by- 
two  matrices  along  the  main  diagonal;  in  the  latter  case  the  covariance  matrix  would  be  a 
general  nonsingular  symmetric  matrix,  possibly  with  relatively  large  off-diagonal  entries. 

In  any  event  the  multivariate  covariance  structure  underlying  the  posit!  >ns  constituting 
the  tracks  must  be  known  and  is  to  be  entered  as  data.  The  scheme  is  based  on  the 
assumption  that  the  degree  of  overhang  for  overlapping  tracks  is  limited  to  one  point  at 
each  end  of  the  tracks;  that  is,  when  two  tracks  start  at  different  times,  the  ^arly  track 
contains  at  most  one  early  point  and  similarly  a track  which  continues  after  the  other 
has  stopped  contains  at  most  one  late  point.  Finally  the  scheme  will  test  track  associa- 
tions for  track  positions  observed  either  at  identical  times  or  at  different  times. 


PROGRAM  STRUCTURE 

The  following  steps  comprise  the  major  activities  of  the  track  association  program. 

Read  Inputs.  For  each  point  in  each  track  the  program  first  obtains  the  time  of 
observation,  the  X coordinate,  and  the  Y coordinate.  In  addition,  it  obtains  the  covari- 
ance matrix  related  to  the  set  of  points  constituting  each  track.  Finally  it  generates  a 
column  vector  TRK  for  each  observed  track,  defined  such  that  for  the  fifth  observed 
track  (K  = 1,  2)  the  /th  observed  (A,  Y)  coordinate  pair  is  defined  as  the  pair 

(TRK  (K,  21  - 1),  TRK  (K,  21)). 

The  times  T used  in  the  program  are  assumed  to  be  defined  in  terms  of  decimal  hours;  some 
changes  will  be  required  whenever  actual  inputs  are  given  in  terms  of  hours,  minutes,  and 
seconds. 

Compute  the  Interpolation  Parameters.  The  fundamental  concept  of  the  program  is 
the  coordinate-by-coordinate  comparison  of  a pair  of  standardized  interpolated  tracks.  The 
program  obtains  these  from  the  track  vectors  TRK(l)  and  TRK{ 2)  by  a linear  interpolation 
scheme  which  produces  a set  of  identical  time  points  and  related  positions  within  each  track, 
The  interpolation  scheme  proceeds  as  described  in  the  main  body  of  this  report  in  conjunc- 
tion with  Fig.  6 and  7.  Next,  if  the  starting  time  for  one  interpolated  track  is  earlier  than 
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that  for  the  other  one,  the  program  deletes  the  earlier  one,  so  that  both  interpolated 
tracks  will  now  start  at  the  same  time;  similarly  any  overhang  at  the  termination  of  the 
tracks  is  reduced  to  the  case  of  equal  track  termination  time.  Thus  both  interpolated 
tracks  will  be  based  on  the  same  set  of  observed  time  points.  The  number  of  such  points 
is  denoted  NPTS;  since  each  observation  point  is  related  to  two  coordinates,  the  majority 
of  the  subsequent  computation  is  based  on  twice  the  value  of  this  parameter. 

Compute  the  Transformation  Matrices  for  the  Original  Tracks.  The  interpolated 
tracks  used  for  the  statistical  test  are  obtained  from  the  original  tracks  by  linear  trans- 
formations. If  the  input  track  is  based  on  M observed  prints  (2 M coordinates),  the  trans- 
formation matrix  Q has  2 (NPTS)  rows  and  2M  columns. 

Compute  the  Interpolated,  Tracks.  For  each  K = 1,  2,  input  track  TRK(K)  is  trans- 
formed into  interpolated  track  TINT(K)  via  the  transformation  matrix  Q,  with  suitable 
dimensions.  Since  the  input  track  is  treated  as  a column  vector,  the  equation  for  this 
transformation  is 


if 

) 

I 

If 

t 

i 


TINT  * Q • TRK. 

Compute  the  Covariance  Matrices  for  the  Interpolated  Tracks.  The  covariance  matrix 
£77  underlying  the  related  interpolated  track  is  given  by 


L -«Z>. 


t where  is  the  covariance  matrix  for  an  originally  obtained  track  and  the  prime  repre- 

i sents  matrix  transpose. 


; Compute  the  Difference  Vector  for  Interpolated  Tracks.  The  hypothesis  that  the 

*.  two  original  tracks  are  from  the  same  target  is  equivalent  to  the  hypothesis  that  the 

? underlying  mean  vectors  are  equal,  or  that  the  difference  of  the  mean  vectors  contains 

y all  zero  coordinates.  This  formulation  carries  over  to  the  interpolated  vectors  as  well, 

| and  at  this  step  the  program  computes  the  value  of  the  difference  TDIFF  of  the  com- 

puted interpolated  tracks,  TDIFF  being  a column  vector. 


$ 

jf 


Compute  the  Matrix  Sum.  In  addition  to  the  values  of  the  difference  vector,  the 
test  requires  knowledge  of  the  covariance  matrix  associated  with  tire  difference,  in  this 
case  the  sum  2 of  the  covariance  matrices  ~LTj  associated  with  the  interpolated  tracks. 

Compute  the  Measure  R 2.  /t2  is  calculated  by  the  formula 


r2  = (TDIFF)' T,-1  (TDIFF), 
where  again  the  prime  indicates  matrix  transpose. 

Compute  Allowable  Upper  Limit  for  R 2.  The  program  next  computes  the  upper 
limit  x2,  above  which  computed  values  of  R 2 will  lead  to  rejection  of  the  hypothesis 
that  the  two  tracks  are  from  the  same  target.  This  limit  is  the  solution  of  expression 
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(4)  of  the  main  body  of  the  report.  These  computations  use  two  approximations  The 
first  one,  based  on  equation  (28)  in  Chapter  17  of  Ref.  10  (Vol.  1),  permits  approxima- 
tion of  the  chi-square  cumulative  distribution  function  (c.d.f.)  by  the  distribution  func- 
tion of  an  appropriately  defined  normal  distribution.  The  second  one,  based  on  expres- 
sion 26.2.23  of  Ref.  16,  uses  the  ratio  of  two  polynomials  to  estimate  percentiles  of 
normal  distributions.  This  method  requires  relatively  little  execution  time  and  produces 
results  which  are  within  a few  percentage  points  of  values  in  the  standard  statistical  tables. 

Make  Decision.  The  program  compares  the  computed  value  of  12 2 with  the  upper 
limit  x2-  II  R2  ^ x2>  the  program  sends  a message  to  merge  the  two  tracks.  Otherwise, 
it  sends  a message  not  to  merge  the  tracks  based  on  the  available  data. 


PROGRAM  LISTING 

Appendix  B is  a listing  of  the  program,  as  written  for  implementation  on  the  CDC 
KRONOS  time-sharing  computer  system.  In  the  KRONOS  version  the  array  sizes  reflect 
the  small  number  of  data  points  used  for  the  test  runs;  for  actual  data  the  array  sizes 
may  have  to  be  increased.  Also,  the  following  items  are  specific  for  KRONOS  and  will 
have  to  be  changed  for  other  modes  of  operation,  such  as  batch  processing  or  implemen- 
tation within  an  automated  correlation  processor: 

• The  array  W1  and  W2; 

• The  CALL  MATIDN  instruction  (line  00160); 

• The  instruction  W2(l,  J)  - S1GM(1,  J ) (line  03320); 

• The  CALL  MATIN V instruction  (line  03660). 


ADDITIONAL  PROGRAMS 

An  extended  version  of  this  program,  which  uses  Kalman  filtering  techniques  to  ob- 
tain smoothed  position  estimates  and  the  related  covariance  matrix,  has  been  developed. 
Information  on  the  extended  program  and  on  other  programs  related  to  the  track  associa- 
tion problem  can  be  obtained  from  the  author. 
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PROGRAM  LISTING 


LIST 

7 V 03/20.  11.28. £4. 
PROGRAM  ASSOC 


k 

'K' 


ft 


I 


| 

§ 


00020  PROGRAM  ASSOC  (OUTPUT.  TAPED 
0 0 04  0 II I MENS  I ON  BUM  <16.16).  M <2>  . HUM  <£)  , 

0 0 06 0 DIMENSION  SI G <2 .16.1 6) . S IGI <16* 1 6) t 
0 0 08  0 Ii  I MENS  I ON  TCOM  <8>  > TDIFF  <16)  > TINT  <2 
0 0 1 0 0 D I MENS  I ON  TRK  <2 , 8 ) , TS  <2 . 8>  . X <2  »4)» 
0 0 1 £ 0 D I MENS  ION  W 1 < 1 6 . 32)  . I..J2  < 16.16? 

00140  ALPHA* 0.05 
0 0 1 6 0 CALL  MAT I DN <W2 .16.1 6> 

00130+ 

002 00+  INPUT  DATA 
0 022 0 READ  <1.5 05 0>  M < 1 > 

0024  0 5050  FORMAT  (13) 

00260  M2=M (1) 

00280  M22=M2+M2 
00300  DO  20  1=1 .M2 

0 032 0 READ (1.5  07  0) TS (1 . I > » X ( 1 . 1) . V ( 1 . I > 
00340  5070  FORMAT (6F8. 4) 

00360  £0  CONTINUE 
00380  DO  30  1=1 . M22 

0 04 0 0 READ (1.5 07  0)  <S6 <1.1. J) . J= I . M22) 

00420  30  CONTINUE 
00440  DO  50  1=2. M22 
00460  11=1-1 
00480  DO  40  J=1.I1 
0 05  0 0 SG  <1.1. J) = JG < 1 . J . I > 

0 0520  40  CONTINUE 
00540  50  CONTINUE 
0 056 0 READ  <1.5 05 0)  M <2> 

00580  M2=M<2) 

00600  M2£=M2+M£ 

00620  DO  60  1=1. M2 

00640  READ  <1 » 5070)  TS  <2. I) . X <2.  I) » V <2, I ) 
00660  60  CONTINUE 
00680  DO  70  1=1. M22 

0 07  0 0 READ  <1.5 07  0)  <S6 <2 . I . J) . J= I . M22) 

00720  70  CONTINUE 

00740  DO  90  I =2 . M22 

00760  11=1-1 

00780  DO  80  . J = 1 . II 

0 0 8 0 0 S G < 2 . I . J ) = ■'  G < £ ; J . I ) 

00820  80  CONTINUE 
00840  90  CONTINUE 
00860  DO  100  K=1 ? £ 

00884  M2=M <K) 

00900  DO  95  1=1 . M2 
0 092 0 TRK  <K . I + 1 - 1 ) =X  <K . I ) 

0 094  0 TRK  <K . I + p = Y <K . I ) 

00960  95  CONTINUE 
00980  100  CONTINUE 


G <2 . 16.  1 6)  . 
S IGM  <16 . 16 
.16).  TM  <2 . 

Y <2 , 4) 
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01000  PRINT  9000 

01080  9000  FORMAT  (30X.8H1.  INPUTS?* 

01040  DO  140  K=  1 > £ 

01060  M2=MOO 
01080  M££=M£+M£ 

01100  PRINT  9010?  K 

01120  9010  FORMAT  9HTRACK  NO . . 1 2 ••'3K . 4HT I ME » 7X . 7HX-C00RD  > 

01 140+5X. 7HY-C00RD) 

01160  DO  130  1=1. M2 

01180  PRINT  9020.  TS(K.I).  X<K,n>  Y(R,I> 

01200  90£0  FORMAT  <5 <F 1 0 , 4 . £X>  ) 

01220  130  CONTINUE 
01240  PRINT  9030.  K 

01260  9030  FORMAT  </•-'.£ OX .2 1 HCOVAR I ANCE  MATRIX  NO.  > I£> 

01280  DO  135  I — 1 » M££ 

0 1 3 0 0 PR I NT  9 04  0 . (SG <K ? I . J) . J= 1 < M22) 

0 1 32 0 9 04  0 FORMAT  <& <F 1 0 . 4 . 1 X) > 

01340  135  CONTINUE 
01360  140  CONTINUE 
0 1 38 0* 

01400*  COMPUTE  THE  INTERPOLATION  PARAMETERS 
01420  CALL  ORDER <TS* M.T.N) 

01440  DO  160  K*1 > 2 
01460  M 1 =M  OO - 1 
01480  NUM <K)  =0 
01500  DO  155  1=1 »N 
01520  DO  150  J=1 . Ml 
01540  IF  ( <T  < I ) . 6E . TS ( K . J)  > . AND . (T  ( I > . LE . T S (K » )+ 1 ) > . 60  TO  145 
01560  GO  TO  150 

0 1 58  0 145  TM  <K . I . J)  = <TS  <K . J+l > -T  >11  > V <TS  <K , J+ 1 > -TS  (K  > J > ':> 

01600  TM  <K.  I > J+l)  = <T  <I>  -TS  <K . J)  > / (TS<K> J+l) -TS  <K»  J)  '> 

0 1 620  NUM  <K> =HUM <K)  + 1 
01640  GO  TO  155 
01660  150  CONTINUE 
01680  155  CONTINUE 
01700  160  CONTINUE 

0 172  0 IF  <TS  (1.1 ) . EQ . TS <2 . 1 > ) GO  TO  32 0 

01740  I F a ( 1) . EG . TS  (1 » 1 > ) GO  TO  200 

01760  1=2 

01780  GO  TO  £40 

0 1 8 0 0 2 0 0 I = 1 

01820  GO  TO  £40 

01840  £40  K-N-l 

01860  NUM  < I > =MUM < I) - 1 

0 188  0 DO  3 0 0 J = 1 « K 

01900  DO  280  1=1.2 

01920  MX=M  (I) 

01940  DO  £60  L=1.MX 
01960  TM  ( I . J. L) =TM (I . J+l . L) 

01980  260  CONTINUE 
02000  £80  CONTINUE 
02 02 0 TCOM  (J)  =T  ( J+ 1 '> 

02040  300  CONTINUE 

02  06 0 NPTS=M  I N 0 (NUM  ( 1 ) . NUM  (£>  ') 

02080  GO  TO  400 

02100  320  NF'T  S=M  I N 0 (NUM  ( 1 > .NUM(£:>) 

02120  DO  360  J=1.NPTS 
0214  0 T COM  < J) =T ( J) 

02160  360  CONTINUE 
02180  400  CONTINUE 
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02200* 

02220*  CONSTRUCTION  OK  THE  MATRICES  Q1?G£  FOP  TRANSFORM INi 
02240*  TRACKS  1 ANH  £> RESPECT  I VELY. 

0226 0 HO  740  K=1.2 
0228  0 M£=MOO 
02300  DO  7£0  I=1,NPTS 
02320  DO  700  -J=l  >M£ 

02340  I NEW= I + 1 - 1 
02380  JNEW=J+ J-l 
02380  CKK> I NEW  j JNEW  > =TM (K « I? J> 

02400  I NEW= I + I 
02420  JNEW=J+J 

02440  O <K> I NEW. JNEW> =TM CK* I » J) 

02460  700  CONTINUE 
02430  720  CONTINUE 
02500  740  CONTINUE 
0252 0* 

02540*  COMPUTING  THE  INTERPOLATED  TRACK  POSITIONS. 

0256 0 NROW=NPT S+NPTS 


0 2580  DO  8 0 0 K = 1 « £ 
02600  NCOL=M  <K> +M  <K> 
02620  DO  780  I=l.NROW 
02640  TINT  <K» I>  = 0. 
02660  DO  760  J=1 . NCOL 


SO 

TINT 

<K, I>=TINT 

00 

7 6 0 

CONTINUE 

£0 

780 

CONTINUE 

40 

8 0 0 

CONTINUE 

60+ 

♦TPk  i'K>  J.)  ) 


02780*  COMPUTATION  OF  S 1 G=C'*SG*  •Q-TRANSPOSE'> 
02800  DO  9 00  K=1  < 2 
0282 0 NCOL=M CK> +M  <K  » 


02840  DO  880  1=1 . NPOW 
02860  DO  860  J = 1 . NPOW 
0 £ 8 8 0 S I G < K j I > J ) = 0 . 

029  0 0 DO  840  N= 1 « NCOL 
0292 0 DUM C I j N> = 0 . 

02940  DO  820  L=t .NCOL 

0296  0 DUM  < I > N>  =DUM  I . N > + (O  f.K . I > L > *SG  '..K . L > N > > 
02980  820  CONTINUE 

0 3 0 0 0 S I G <K . I > J > =S  1 6 <K » I . J>  + < DUM  < I . N.)  *0  Of . J . N > > 

03020  840  CONTINUE 

03040  860  CONTINUE 

03060  380  CONTINUE 

03 08 0 900  CONTINUE 

03100  NCOL=NROW 


v,  . 
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031 20*  •*• 

03140*  COMPUTE  TD I FF =T I NT  1 -T I NT2  > DIFFERENCE  VECTOR 

03160  DO  320  l=l.NROW 

03180  TD IFF ‘ I) =T I NT <1. 1 ) -TINT <2.  II 

03200  920  CONTINUE 

03220* 

03240*  COMPUTE  SH3M=SIGl+5 162.  MATRIX  SUM 
03260  DO  960  1=1. HROW 
03280  DO  940  J=1 . HCOL 
0 33  0 0 1 1 bM  ( 1 . . J > —13 1 b *i.  1 . I . J !'  +3 I b <2  > I . J.) 

0 332  0 we  < I , J)  =3 1 6M  (.1  > Jj 
03340  940  CONTINUE 
03360  960  CONTINUE 
03380  PRINT  9050 

03400  9050  FORMAT  <■•••  -•'>  20X«  20H2.  TRANSFORMED  VALUES  ^ 1 9X . 7HTRACK  Is 
0342  0+ 1 ?X - 7HTRACK  2/3X.  4HT I ME » 2X » 2 v.5X . 7HX-C00RD . 5X  * 7HY-C  C10PD)  ) 

03440  DO  965  1=1 .HP T 3 

03460  PRINT  9020.  TCOM  <1  ) . TINT  Cl . 2*1 -1>  .TINT  <1.2*0  .TINT  <8»g*I-l)  > 
0348  0+T I NT  (2.2*0 
03500  965  CONTINUE 
03520  PRINT  9060 

03540  9060  FORMAT  <••  -•.  15X. 38HSUM  OF  TRANSFORMED  COVARIANCE  MATRICES) 
3560  DO  970  I=l.NROW 

0 358 0 PRINT  9 04 0 > C S I GM < I . J) . J= 1 , NPDM) 

03600  970  CONTINUE 
0362 0* 

03640*  COMPUTE  P2=  (TD IFF -TRANSPOSE) ♦ c S I GM- INVERSE) ♦ <TD IFF) 

03660  CALL  MATINV <W2. SIGI . 16.32.W1) 

>1368  0 R 2 = 0 » 

0 3700  DO  1 000  J=  1 . NF'QW 
03720  DIJM  (1 . J)  =0. 

03740  DO  980  I=l.NROW 

0376  0 DUN  < 1 . J)  =DUN  < 1 . J">  + (TD  I FF  < I » *S  16 1 < I . J. ) 

03780  930  CONTINUE 

038.00  R£=R2+  (HUM  1 1 > J)  *TDIFF  < J)  > 

03820  1000  CONTINUE 
03340  PRINT  8000.  R2 

03860  8000  FORMAT  W.20X.  1 7H3.  COMPUTED  VALUES. , 3HP2=.  F12.5) 
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03880* 

03900*  COMPUTE  THE  <1 -ALPHA? -TH  PERCENTILE  POINT  FOR.  THE 

03920*  CHI-SQUARE  DISTRIBUTION  WITH  2*NPTS  DEGREES  OF  FREEDOM 

03940  XNPT=2*NPTS 

03960  PROB=l.- ALPHA 

03980  IFCPROB.LE. 0.5>  GO  TO  1020 

04000  P=1 . -PROB 

04020  GO  TO  1040 

04040  1020  P=PROB 

04060  1040  TE=SQRT CRLOG <1 . / <P*P? > > 

04080  C 0=2. 515517 
04 100  C 1 = 0 . 8 02853 
04120  02=0.010328 
04140  01=1.432788 
04160  D2= 0.189269 
0 4180  D 3 = 0 . 0 0 1308 

04200  U 1 = 1 . + <D1*TE>  + <D2*TE*TE)  + <D3*TE*TE*TE? 

04220  U2=C0+ <C1*TE>  + <C2*TE*TE> 

04240  XP=TE~ <U2/U1> 

0426 0 XF-XP*  <PROB-0 . 5> CABS  CPROB- 0 . 5)  > 

04280  A=  CXP*SORT  <2.  <9.  *XHPT>  » > +1 . 

043  0 0 CH I SG=XNPT*  < <A-  <2 . < 9 . *XNPT?  > ')  ♦ *3> 

04320  PRINT  8020 . ALPHA,  XNF'T . CHI  SO 

04340  8020  FORMAT  C6HAIPHA=>  F6.  3-  2X,  15HDEG.  OF  FREEDOM. F4. 0. 
0436  0+2X » 6HCH I SQ= . F 1 2 . 4 :• 

04380* 

04400*  DECISION  POINT 

0442  0 I F CR2-CH I SQ>  1 06  0 « 1060*  1 08 0 

04440  1060  PRINT  8040 

04460  804  0 FORMAT  <1 2HMERGE  TRACKS'* 

04480  GO  TO  2000 
04500  1080  PRINT  8060 

04520  8060  FORMAT  C32HDO  NOT  MERGE  TRACKS  AT  THIS  T IME* 

(14540  GO  TO  2000 
04560  2000  CONTINUE 
04580  STOP 
04600  END 


i 

l 

i 


I 


i 

I 


■j 
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04620 
04640 
04660 
04680 
047  0 0 
04720 
04740 
04760 
04780 


SUBROUTINE  ORDER <G > MM* H* N> 

D I MENS  I ON  0 <2  * 1 5)  > MM  <2')  * H < 15  > > XX  ( 1 5) 
MMI=MMa> 

DO  20  I = 1 * MMI 
XX  < I ) =©  ( 1 > I 
£0  CONTINUE 
NMM=MM  < 1> +MM C2> 

MH2=MM C2> 

DO  40  1 = 1*  MM2 


04800  J=MM a> + 1 
04820  XX  CJ ' =Q  <2 * I > 
04Q40  40  CONTINUE 
04860  DO  80  _l=l * MMM 


04880  H'-J'=XX(J> 

04900  DO  60  I=J*MMM 

04920  IF  <XX<1> .GE.HfJV*  60  TO  60 

0494  0 Z=H  (J.) 

0496  0 H < J > =XX  < I > 

04980  XX<I>a2 

05000  60  CONTINUE 

05020  80  CONTINUE 

05040  N=NMM 

05060  DO  120  1=1 -MMM 

05080  IF  <Ha>  .NE.H.\I  + r>  ' 60  TO  120 

05100  11=1+1 

05 1 £ 0 DO  1 0 0 J=I 1 ’ MMM 

0514  0 H >’  J)  =H  < J+l  > 

05 1 6 0 1 0 0 CDNT I NUE 


05180  N=N-1 
05200  MMM*MMM- 1 
05220  120  CONTINUE 
0524  0 MMM=MM  < 1 > +MM  <.£> 
05260  RETURN 
05280  END 
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UNITED  STATES  GOVERNMENT 

memorandum 

date:  5 January  1996 


REPLY  TO 

ATTN  OF:  Code  8150 

SUBJECT:  REVIEW  OF  NRL  REPORT  7989 


Code  1221 


1.  It  was  requested  that  I review  NRL  Report  7989,  " Track- to-Track  Association  in  Ocean 
Surveillance",  by  Howard  Weiner  for  a possible  change  in  distribution  statement.  I have 
subsequently  reviewed  the  copy  provided  and  can  see  no  reason  for  restricting  access  to  this 
information. 


2.  It  is  recommended  that  the  distribution  statement  on  this  subject  report  be  changed  to 
unlimited  public  release. 


R.L.  Beard 

Head,  Space  Applications  Branch 
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